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Abstract 

We study a linear threshold agent-based model (ABM) for the spread of political revolutions on 
social networks using empirical network data. We propose new techniques for building a hierarchy of 
simplified ordinary differential equation (ODE) based models that aim to capture essential features of 
the ABM, including effects of the actual networks, and give insight in the parameter regime transitions 
of the ABM. We relate the ABM and the hierarchy of models to a population-level compartmental 
ODE model that we proposed previously for the spread of political revolutions pQ, which is shown to be 
mathematically consistent with the proposed ABM and provides a way to analyze the global behaviour of 
the ABM. This consistency with the linear threshold ABM also provides further justification a posteriori 
for the compartmental model of pQ. Extending concepts from epidemiological modelling, we define a 
basic reproduction number Ro for the linear threshold ABM and apply it to predict ABM behaviour 
on empirical networks. In small-scale numerical tests we investigate experimentally the differences in 
spreading behaviour that occur under the linear threshold ABM model when applied to some empirical 
online and offline social networks, searching for quantitative evidence that political revolutions may be 
facilitated by the modern online social networks of social media. 


1 Introduction 

In this paper, we propose a linear threshold agent-based model (ABM) for the spread of political revolutions 
on social networks in dictatorial regimes. In our ABM, nodes of the network can be in two states: active, i.e., 
participating in the revolution, or inactive. Transitions from the inactive to the active state are governed by 
a growth process that uses a traditional linear threshold mechanism: an individual v may change from the 
inactive to the active state if the fraction of neighbours in the active state exceeds the linear threshold 6. 

Linear threshold models have been studied before in the context of influence maximization [2] , enforcement 
of unpopular norms (3j, and general “complex contagion” processes PI The development and analysis of 
linear threshold processes has been highly influenced by more established epidemiology mm and rumour 
spreading nm nu models, with one key difference: whereas a contact with a single “infected” individual 
is enough to spread a contagion in most epidemiology or rumour spreading applications, i.e. in “simple 
contagion” processes, multiple contacts with infected individuals are generally necessary to spread a contagion 
in a linear threshold process. Indeed, epidemiological and rumour spreading models can be considered to be 
a special case of linear threshold models where the threshold parameter is chosen sufficiently small such that 
one infected neighbour enables propagation. Most of the linear threshold models that have been studied in 
the context of social spreading or contagion processes consider evolution in discrete time. In 

contrast, our ABM is chosen to evolve in continuous time since that facilitates the comparison of the ABM 
with simplified population-level models, which is a major goal of this paper. 

The main justification for our linear threshold modelling approach in the context of political revolutions 
is that the decision to join a political revolution can be assumed to be a collective action problem US]: if 
individuals act unilaterally they are subject to retaliation by the regime, whereas if they act collectively then 
the regime loses the ability to punish due to a lack of resources. The linear threshold modelling approach is 
consistent with the collective action principle, since an individual transitions from an inactive to an active 
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state only after the number of neighbours in the active state has been observed to reach the critical fraction 
9 , and the individual deems the revolution of sufficient size to consider participation. In order to incorporate 
communication network structure explicitly, we represent interactions between individuals v by edges e in an 
undirected graph G = G(V,E). This contrasts with previous attempts to model political revolutions using 
linear threshold models that assume homogeneous mixing mug, that assume that G — G(V, E) is the 
complete graph. 

The goals and contributions of this paper are as follows. We start by presenting the linear threshold ABM 
for the spread of political revolutions on social networks, and then derive a hierarchy of simplified ordinary 
differential-equation (ODE) models of varying degree of sophistication that characterize the solutions of the 
linear threshold ABM, see Table |T| Although the linear threshold ABM allows us to explicitly model the 
effects of the communication network on the dynamics of the political revolution it models, the inherent 
complexity of modelling the state of all individuals simultaneously makes this approach difficult to analyze 
and expensive to simulate. For these reasons there is significant interest in simplified aggregate or population- 
level dynamical models, which can provide more cost-effective dynamical simulations, and can give insight 
in the qualitative dynamics of the ABM and its parameter regimes, through dynamical analysis of the 
simplified model. It is a significant challenge to incorporate actual network structure in the ODE models, 
and this is often important since details of the network structure may determine the qualitative behaviour 
of ABM solutions as the ABM parameters are changed. In this paper, we present two new effective ways to 
incorporate network structure into the one-compartmental ODE that approximates the dynamical evolution 
of the expected fraction of the population that participates in the revolution in the linear threshold ABM 
model. These approaches make use of the degree distribution of the graph or samplings of the network, and 
we call the resulting population-level ODE models the binomial visibility function (BVF) and the empirical 
visibility function (EVF) models. The EVF model is significantly less expensive than the BVF, but we 
demonstrate its equivalence with the BVF in the limit of large network and sample size. Extensive numerical 
tests show that these ODE models perform remarkably well to predict qualitative behaviour of the ABM; 
in fact, they often perform as well as or better than a previously specified higher-order model that has 
many compartments and is much more expensive computationally |l4j . and that we refer to as the degree 
approximation (DA) model. 

Table 1: Summary of hierarchy of models from highest (network-level) to lowest (population-level) complex¬ 
ity. 


Model Abbreviation Section Equation(s) Reference(s) 

Linear Threshold Agent-Based Model ABM 

Degree Approximation Model DA 

Binomial Visibility Function Model BVF 

Empirical Visibility Function Model EVF 

Step Visibility Function Model SVF 
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Similar to the population-level one-compartmental ODE model that we proposed in [lj for the spread of 
a political revolution, the linear threshold ABM proposed in this paper applies to dictatorial regimes that 
employ censorship and police repression. We make the simplifying assumption that the population is uniform 
in its dislike of the regime (such that all nodes can potentially become active in the revolution), and we assume 
that the regime censors the communication between individuals except for the links in the social network. 
A decay process models police repression that is limited by the regime’s finite police capacity. The model 
features parameters that describe timescales of growth and decay, a linear threshold parameter for the growth 
process, and a parameter for the finite police capacity. The dynamics of our ABM is governed by stochastic 
transitions in continuous time, because this allows us to relate the ABM mathematically to simplifying time- 
dependent ODE models in limits of expectation over the aggregate population. These simplified ODE-based 
models aim to capture essential features of the linear threshold ABM at the population level, including effects 
of the actual networks, and give insight in the parameter regime transitions of the ABM. 

We also relate the ABM to the population-level one-compartmental ODE model for the spread of political 
revolutions that was proposed in [1] and that indirectly takes network structure into account by a single 
parameter that describes from which fraction of participants the growth term is switched on discontinuously. 
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As such, it is a precursor of the BVF/EVF models introduced in this paper. In pLj this single parameter 
was called the visibility parameter of the revolution, with the interpretation that people will only join the 
revolution when the revolution is of sufficiently large size to be visible to the population, while repressive 
regimes attempt to make unrest invisible through censorship. We call the model from |Tj the step-visibility 
function (SVF) model in this paper. This ODE model was justified in JT] by a simplified network model 
with significant further assumptions; in this paper, we show that our ABM is mathematically consistent 
with the model from [T| and we test the assumptions from |pQ with the BVF/EVF models and the ABM 
model on real networks, corroborating a posteriori the assumptions from pQ, and showing that the model 
from [T] reproduces the qualitative behaviour of the ABM and the spread of a revolution under a linear 
threshold model. This provides further justification for the ODE model of pQ, and at the same time provides 
inexpensive ways to predict the parameter regime behaviour of the ABM presented in this paper using 
simple compartmental models like the SVF and the new BVF/EVF models. We emphasize that the broad 
applicability of the linear threshold process implies that the framework, methods, and results presented in 
this paper are relevant to the study of many related social spreading processes. 

The final contribution of this paper relates to an interesting application of our linear threshold ABM. We 
investigate experimentally the differences in spreading behaviour that occur under the linear threshold model 
when applied to some empirical online and offline social networks, searching for quantitative evidence that 
political revolutions may be facilitated by the online social networks of social media. Indeed, it is often 
assumed that the connectivity of modern online social networks has greatly facilitated the spread of political 
revolutions in the past decade, e.g., in the Arab Spring revolutions of 2011 ns ceei nzi m m, while the 
traditional offline social networks of the pre-Internet era (using in-person physical contact, or mail or phone 
interaction for safe communication) had a different connectivity structure that was often severely restricted 
by regime censorship. As a starting point to investigate differences in propagation properties that may arise 
between online and offline social networks within the linear threshold propagation model, we investigate 
propagation on two empirical networks of modest size using our linear threshold ABM: we consider a small 
Facebook network m as a representative of an online social network, and a small physical contact network 
between individuals |20j as a representative of an offline social network. These two networks are also used 
for all validation and comparison simulations for the ABM, DA, BVF/EVF and SVF models throughout the 
paper. We also extend the concept of basic reproduction number R 0 from epidemiological modelling [9] to 
characterize the potential of networks to spread the revolution. The basic reproduction number for our ABM 
with a linear threshold process is easy to compute, and we show that it gives useful predictions. Our results 
indicate that the offline social network indeed is less conducive to spreading the revolution than the online 
social network: it has a smaller basic reproduction number, and in simulations a larger initial population 
of revolutionaries is required to spread the revolution. This provides some initial quantitative evidence that 
the spread of revolutions under a linear threshold process may occur more easily on modern online social 
networks than on traditional offline networks, but we also comment on limitations in our approach and further 
investigations that are required to address this intriguing but complex question more comprehensively. 



(a) Physical contact network ( Gp ) (b) Facebook social network (Gp) 


Figure 1: Visualization of physical contact (Gp) and Facebook (Gp) social networks using the Gursoy-Atun 
algorithm m as implemented in the Matlab BGL package. 
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The rest of this paper is structured as follows. Section [2] describes the Facebook and physical contact 
empirical networks on which we test our ABM and approximations. Section [3] specifies the linear threshold 
ABM in detail and briefly reviews the SVF model presented in [T]. Section [3] then establishes the consistency 
between these two approaches through mathematical derivation and numerical simulation on the networks 
described in Section [2] i.e., we show that the SVF model can be considered to approximate the aggregate 
behaviour of the ABM. Section [4] describes the BVF and EVF methods to improve the SVF model by 
changing the functional form of the visibility function, explicitly incorporating the structure of the network 
without increasing the dimensionality of the SVF model. Section [5] confirms the usefulness of the BVF and 
EVF models by showing that they perform no worse than the higher order degree approximation model 
of dl], while being much less expensive to evaluate. Section [h] extends the concept of basic reproduction 
number from epidemiology [9J to the linear threshold process specified in Section [3] and illustrates how it can 
be interpreted in terms of the BVF/EVF. Finally, Section [7] explores how our methods can be used to study 
the differences in spreading behaviour of political revolutions on online and offline social networks under the 
linear threshold mechanism. 


2 Network Data 

In this section we describe the two empirical social networks that will be used in this paper to validate and 
compare the ABM, DA, BVF/EVF and SVF models for the spread of political revolutions on social networks 
under the linear threshold model. 

The first network we consider is the physical contact network presented in |20j . It was constructed by 
distributing wireless sensors to students, teachers, and staff at a U.S. high school during a one day period. 
When two wireless sensors are in proximity of one another, i.e. when they are less than approximately 3m 
apart, they register an interaction with a temporal resolution of 20s. Therefore, the communication network 
we extract from this data is referred to as the physical contact network Gp. The second network we consider 
is the Facebook subnetwork presented in [IS]. It was constructed by combining the ego-network ^ of ten 
individuals and then taking the largest connected component of the resulting network. We refer to this 
network as the Facebook social network Gf- For details on the network extraction protocol for Gp we refer 
to Appendix[A] To facilitate the comparison of experimental results, the minimal contact duration to register 
an edge in Gp is chosen such that Gp and Gp have approximately the same average degree, see Appendix 

El 

The physical contact network Gp has N = 776 nodes with sparsitj^S 1 2 3 = 0.06, and is visualized in Fig. 

The Facebook network Gp has N = 3,963 nodes with sparsity S = 0.01, and is visualized in Fig. |lb| The 
cumulative degree distribution^] are displayed in Fig. [ 2 ] The distributions of local clustering coefficients are 
displayed in Fig. [3] 

It is interesting to consider how the physical contact and Facebook networks differ on some common network 
measures, see Table [2} Comparing the physical contact and Facebook social networks along these admittedly 
limited measures nevertheless highlights substantial structural differences between these networks. The 
physical contact network appears more homogeneous than the Facebook social network, in the sense that 
the Facebook social network can be visually grouped into distinct communities, whereas the physical contact 
network cannot (Fig. [I]). The network measures presented in Figs. [ 2 J 3 ] and Tableware generally supportive of 
this observation. For example, the cumulative degree distributions displayed in Fig. [2]show that the physical 
contact network has an approximately normal distribution with a relatively small standard deviation, and 
hence, a thin tail. In contrast, the Facebook social network has a much broader distribution (possibly scale- 
free or exponential, although this cannot be determined conclusively with such a small network size) with a 

1 For a network G = G(V,E) the ego-network of an individual v S V is a subnetwork of the overall network composed of 
individual v (called the ego-node), the neighbours of individual v, and all of the connections between these individuals. 

2 The Sparsity of a network is defined as the fraction of possible edges that are present in the network. 

3 If f>k is the fraction of nodes with degree k (i.e., the degree distribution) of the graph, then the cumulative degree distribution 
is the function F(k) = )T/’ = 1 Pj- The cumulative degree distribution is displayed in place of the degree distribution because the 
degree distribution is subject to significantly more noise than the cumulative degree distribution. This is a common approach 
for empirical networks. 


la 


4 





Figure 2: Cumulative degree distribution of physical contact (dashed) and Facebook (solid) social networks. 




Local Clust. Coeff. Local Clust. Coeff. 

(a) Physical contact network (Gp) (b) Facebook social network ( Gf) 

Figure 3: Distribution of local clustering coefficient for physical contact (Gp) and Facebook social (Gp) 
networks. 


relatively large standard deviation, and hence, a fat tail. Similarly, it is known that online social networks 
often have larger clustering coefficients than offline social networks. In Section [7] we use the Facebook and 
physical contact networks as representatives of online and offline social networks, respectively, to investigate 
differences in propagation properties that may arise between online and offline social networks within the 
linear threshold propagation model. 
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Table 2: Summary of common network measures for physical contact and Facebook social networks. 


Quantity 

Physical Contact Network 

Facebook Social Network 

Number of Nodes ( N ) 

776 

3,963 

Average Degree (pk ± Ofc) 

44.8 ± 18.2 

44.5 ±52.4 

Minimum/Maximum Degree (fc m inA'max) 

1/109 

2/1034 

Sparsity ( S ) 

0.06 

0.01 

Diameter ( D ) 

8 

8 

Average Path Length(/i; ± 07) 

2.37 ±0.82 

3.77 ± 1.29 

Acerage Local Clust. Coeff. (p c ± <r c ) 

0.29 ±0.12 

0.62 ±0.20 


3 Agent-Based Model 

3.1 Specification of the Agent-Based Model 

Consider a population of individuals who are represented by nodes V = and whose interactions are 

represented by edges E = {e,}-^ 1 , of the graph G = G(V,E) with degree distribution pk- An individual 
v £ V at time t can be in one of two states, i.e. s v (t) = 0 (inactive in the revolution) or s v (t) = 1 (active 
in the revolution). We assume that the network G is static so that the dynamics of the ABM can be fully 
specified by providing rules for the transition of individuals from an inactive state to an active state and vice 
versa, together with parameter values and an initial condition. For the moment, we set aside the issue of 
the choice of parameters and initial conditions and restrict the remainder of this section to the specification 
and justification of the transition rules. 


3.1.1 Growth Process: Inactive ( s v = 0) to Active ( s v = 1) 


To specify the growth process of the ABM, we use the standard linear threshold model that has been 
used before in the context of social spreading or contagion processes such as opinion formation, technology 
adoption, marketing, rioting, and political movements mmmmm- As mentioned in Section[l] the choice 
of individuals to join a revolution is a collective action problem [13] : individuals are averse to unilateral action 
against the regime for fear of severe retaliation, but are willing to take action collectively in the belief that 
the regime will lack sufficient resources to punish the entire collective. Thus, an individual v £ V will decide 
to join the revolution if he or she believes that it has grown sufficiently large to reduce the risk of retaliation 
from the regime to an acceptable level. It is reasonable, therefore, to assume that individuals will decide to 
join the revolution only after a large enough fraction of their neighbours in the social network have done so. 
This behaviour is captured by the following transition rule: if individual v £ V has k neighbours {wj}j =1) 
if s v (t) = 0, and if 

k 

y ^s Wi (t) > e v k, (i) 

J=1 


i.e., if v is inactive at time t and has at least a fraction 9 V of its neighbours that are active, then node v 
transitions from state 0 to state 1 at time t' = i±£i, where > 0 is the first arrival time of a Poisson process 
with rate C\. We say that nodes that satisfy Equation 0 can “see” the revolution, i.e. the revolution is 
visible to them. Alternatively, we say that these nodes are “considering joining the revolution”. Parameter 
ci determines the timescale of the growth process. 


While most of the linear threshold models that have been studied in the context of social spreading or 
contagion processes ® (a no o ni consider evolution in discrete time by choosing £i constant, we have 
specified so that Equation (fT|) evolves in continuous time. Specifying in this way is often considered 
in biology and epidemiology j7J [5] , since it (a) facilitates the comparison of the ABM with population- 
level models (as the SVF model described in Section 3.3), and (b) eliminates the problem of choosing a 
suitable time-step for iterating the discrete-time process. Furthermore, we note that by choosing to be 
the first arrival time of a Poisson process, we are assuming that the decision making process is a Markov, 
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or memoryless, process. Specifically, conditioned on the state of the system at time t, we assume that the 
likelihood of an individual joining the revolution in the time interval [t, t + At] is independent of the state 
of the system at any time t < t. In words, the future is independent of the past, given the present. 

3.1.2 Decay Process: Active (s v = 1) to Inactive (s„ = 0) 

We assume that once an individual v £ V has become active the regime will attempt to arrest or disperse 
him or her, thus returning v to an inactive state. As in the previous section, we assume that this is a 
memoryless process. We further assume that the regime can only arrest or disperse protesters so long as 
the total fraction of active protesters remains less than the regime’s finite police capacity /3 £ (0,1). The 
transition rule can then be characterized by the following: if s v (t) = 1 and 

s w{P) ^ Pi ( 2 ) 

wGV 

i.e., if v is active at time t and the fraction of active individuals is less than the regime’s police capacity /?, 
then the node v transitions from state 1 to state 0 at time t' = t + £ 2 , where £2 > 0 is the first arrival time 
of a Poisson process with rate C 2 , which determines the timescale of the decay process. This fully specifies 
the evolution of the ABM, which we simulate using Gillespie’s algorithm, see Appendix |D| 

3.2 Brief Review of the SVF Population-Level Compartmental Model of [T] 

0 p 1-a 1 

II f-«- fr <-*- 

0 c* 1-a p 1 

1110 {—<-HH—«—I-»- } 

0 1-a c* p 1 

llle i -«- ](> • <)[ - * -* 

0 1-a p c* 1 

mi i—* -3(— * —l-H—*- } 

(b) Equilibria, stability and basins of attraction on the 
r-axis (r £ [0,1]) for parameters a, /?, ci, and C 2 of 
the SVF model lying in Regions I, 1110, llle and III1. 
Closed (open) circles represent locally asymptotically 
stable (unstable) equilibria. Left (right) arrows indicate 
regions where r < 0 (r > 0). 

Figure 4: Summary of the potential dynamical regimes of the SVF model of pQ, given in Equation (|3|, with 
parameters a, /3, c\, and C 2 . 

Before deriving a population-level ODE for the linear threshold ABM and establishing the consistency of the 
SVF model from [T] with the ABM in the following section, it is useful to first review some basic properties 
of the SVF model. For example, we will show later that the parameter regimes of the SVF are predictive of 
the parameter regime transitions in the ABM. The SVF model is given by the equation 

r = ci (1 - r) v s (r;a) — c 2 r p(r;/3), (3) 

v ^ ^ v ^ ^ 

g(r) d(r) 

with parameters a,/3 £ (0,1) and Ci,C 2 > 0. The growth and decay of the fraction of protesters r(t) are 
modelled by the growth and decay functions g , d : [0,1] —> R + , respectively. Given their roles in determining 
the rate of growth and decay, the parameters C\ and C 2 are called the protesters’ enthusiasm and police 


CCl 



c : /{c : +c 2 ) 


a 


(a) Division of a — /3 parameter space of the SVF model 
into Regions I, II, Hie, 1110, and III1. 
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efficiency parameters in [Tj. The growth and decay in the number of protesters are modulated by the 
visibility and policing terms that are taken to be step functions v s (v,a) = I{t.>i_ q } and p(r; /3) = H{ r - </ 3 }? 
where the indicator function is defined as 

j _ f 1 if X is true 
Vd | 0 otherwise 

These functions are parametrized by visibility and police capacity parameters a and (3, respectively. In words, 
the fraction of protesters can grow only when the protest movement is sufficiently large to be “visible” to 
the general population, i.e. when protests are large enough that their existence cannot be masked by the 
regime through censorship. Conversely, the fraction of protesters can decay only so long as the fraction of 
protesters is less than the regime’s finite capacity to disperse or arrest protesters. This is consistent with 
the idea that individuals who consider joining the revolution are faced with a collective action problem Il3l : 
if they choose to act unilaterally then they face severe retaliation from the regime, so they wait to act until 
the revolution becomes visible and if they act collectively then the regime loses the ability to punish due to 
a lack of resources. 

Depending on the values of the parameters a, /3, ci, and C 2 , the dynamics of the SVF fall into one of four 
main regimes, see Fig. [4] The equilibrium of total state control r = 0 and the equilibrium of the realized 
revolution r = 1 are present in all parameter regimes and are always locally asymptotically stable. Different 
parameter regimes are distinguished by the behaviour of Equation ([3]) when r lies between 1 — a and /?. 
Region II (/3 < 1 — a) is characterized by an interval of stable equilibria r £ (/3,1 — a) bounded by two 
unstable equilibria r € {/3,1 —a}. Since parameters in Region II are characterized by a relatively weak regime 
(low /3) and relatively small visibility (low a), Region II is interpreted in [_T| as the failed state parameter 
region. In contrast, Regions IIIO, Hie, and III1 (1 — a < (3) are characterized by a relatively strong regime 
(large (3) and relatively large visibility (large a) and differ in the value of the quantity c* = C\/[c\ + C 2 ), i.e. 
in the ratio of protester’s enthusiam to the sum of protester’s enthusiasm and police efficiency. In Region 

1110, i.e. when c* < 1 — a < (3, the interval (1 — a,/3) lies in the basin of attraction of the equilibrium of 
total state control r = 0. Analogously, in Region III1, i.e. when 1 — a < /3 < c*, the interval (1 — a , (3) lies in 
the basin of attraction of the equilibrium of the realized revolution r = 1. Finally, in Region Hie, i.e. when 
1 — a < c* < /?, the interval (1 — a,/?) becomes the basin of attraction of a new equilibrium r = c*, which 
is called the equilibrium of civil unrest in [Ij. Because of the contribution of the interval (1 — a,/3) to the 
stability of the equilibrium of total state control, of the realized revolution, or of civil unrest, Regions IIIO, 

111 1 , and Hie are called the stable police state , unstable police state , and meta-stable police state , respectively. 
For a full description, interpretation, and analysis of the SVF model we refer the reader to [F. 


3.3 Population-Level ODE for ABM and Relation to SVF Model 

3.3.1 Mathematical Analysis 


To show consistency of the ABM proposed in Section 3.1 with the SVF 
r a (t) = r a {t\t 0 ) to be the fraction of nodes (in the ABM model) that are 
expected fraction of active nodes at time t conditioned on information at 


model of [I] we begin by defining 
expected to be in state 1 , i.e. the 
time to: 


r a{t) = -^X^EK(i)|t 0 ]. 

i= 1 

The fraction of nodes that are expected to be in state 0 (inactive) at time t , conditioned on information at 
time to, is then 1 — r a (t\to). We now write the change in r a from time to to time t as 


Ar a (t) = r a (t) - r a {t 0 ) = g(t\t 0 ) - d{t\t 0 ), (4) 

where At = t — to, and the expected growth and decay of the fraction of active nodes are modelled by the 
non-negative growth and decay functions g(t) = g(t\to), and d(t) = d(t\to), respectively. In order to obtain 
the one-compartment model we will need to approximate the quantities g(t ) and d(t) in terms of r a (to), 
and /3. 



For notational convenience we begin by considering the case where Wv £ V : 9 V = 9. The pool of individuals 
that can go from active to inactive at time to is r a (to). Active nodes become inactive at the first arrival time 
of a Poisson process with rate C 2 , provided that the fraction of active nodes does not exceed the regime’s 
police capacity, i.e. r a (to) < /3. Thus, 

d(t\t 0 ) = [C 2 r a (t 0 ) At + o(At)} p(r a ; fj) = c 2 r a (t 0 ) p(r a (t 0 ); P) At + o(At), (5) 

where we take p(r; j3) = H-[r-c/ 3 > , as in Equation 0. Next, let 


v{r a {to)\6) 


1 

N 


Ato)>Bk v } ’ 

.vev 


( 6 ) 


where denotes the neighbours of individual v £ V. In words, we let v{r a (to)',9) be the fraction of 

the total population at time to that can see the revolution. We call v{r a (to)\ 9) the visibility function of the 
linear threshold ABM. In the chosen notation for v it is emphasized that v depends on r a (to), i.e. on the 
fraction of nodes that are active at time t 0 . At time t = to, the pool of individuals that can go from inactive 
to active, i.e. that are considering joining the revolution, is therefore approximately (1 — r a (to)) v{r a (to)~, 9). 
Since inactive nodes that can see the revolution become active at the first arrival time of a Poisson process 
with rate c\, we have 

g(t\t 0 ) = ci (1 - r a (t 0 )) v{r a {to)',9) At + o(At). (7) 

Combining Equations 0-0 and dividing by At gives 
At 

= a (1 - r a (t 0 )) v{r a (to)\9) - c 2 r a (t 0 ) p{r a (t 0 ); /?) + o(l). 

Taking the limit as At —> 0 yields 


dr 


a 


dt 


ci (1 - r a ) v(r a \9) - c 2 r a p(r a ; f3). 


( 8 ) 


This establishes a general population-level ODE model that approximates the linear threshold ABM, with 
visibility function v{r a \ 9). To close the model, the visibility function has to be specified, guided by Equation 
0. Observe that if we substitute v(r a ;9) by the step function v s (r a ; a(9)), then this exactly recovers the 
SVF model of [T|, see Equation 0 - So, provided that the step function v s (r;a(9)) is a suitable choice for 
v(r\ 9), which we argued in the appendix of ][T] based on considering a simplified averaged threshold process 
that ignores correlations, this confirms the consistency of the ABM and the SVF. The crux of the argument 
used in pQ, which we discuss further in Sec. 4.1 is that, for real networks u(r;9) can be expected to be a 
steep sigmoidal function of r, which can be approximated by a step function. In this paper, we test these 
assumptions explicitly on the empirical physical contact and Facebook social networks discussed above using 
ABM simulations, and we derive new improved ODE models that form a more accurate representation of 
v(r\ 9) than a step function, taking into account the actual network structure of the empirical networks. 
The step function assumption was chosen in pQ because it is useful when little information is known about 
the underlying network structure. For example, without having to specify an underlying communication 
network, the SVF model provides a potential mechanism by which networks with increased visibility a are 
more susceptible to political revolution. 


The following section explores the relationship between the ABM and the SVF model numerically for 
empirical networks by using a fitting procedure to attempt to estimate the network-dependent a = a(9 ) for 
which v s (r ; a(9)) best approximates u(r; 9). Given this estimate for a = a(9), we compare ABM simulations 
to solutions of Equation 0 in order to verify the consistency of the ABM and the SVF numerically and to 
investigate the ABM behaviour in the parameter regimes of the SVF model as described in Fig. [1J Section 
[4] then expands on this by proposing and analyzing a more sophisticated approximation for the visibility 
function v(r~, 9). We conclude this section by noting that the above derivation can easily be modified for the 
case where the 9 V are not uniform. In this case the growth and decay functions g and d, respectively, can be 
computed by averaging over the distribution of the 9 V ’ s. 
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(b) Physical contact network Gp (a = 0.834) 

Figure 5: Parameter Region 1110. (Left) Time traces of ABM simulations with optimally fitted a (grey) and 
solution to Equation <[3j) (black). (Right) L\ error between ABM realizations and the solution to Equation ([3]) 
versus a. The optimal a is found by fitting Equation ([3| (the SVF model) to rep = 100 ABM simulations 
with parameters 9 = 0.15, p = 0.3, Ci = 1, C2 = 9, and ro = 0.25 (Region 1110, c* < 1 — a < P). 


3.3.2 Numerical Simulation 

In this section we simulate the ABM on the physical contact and Facebook networks for a fixed set of 
model parameters and find the optimal a for Equation ([3]) that minimizes the difference between the ABM 
realizations and the output of Equation ([3]). The closeness of the output of Equation ^ and the ABM 
realizations for broad ranges of parameter values demonstrates that the SVF model does indeed approximate 
the ABM, at least for non-pathological parameter choices. We begin by fixing the initial condition r a (0) = 
and the parameters 9 V = 9, ft, c i, and C 2 - To compare realizations of the ABM with the solution to 
Equation © we use Gillespie’s Algorithm (summarized in Appendix |d| to simulate rep = 100 realizations 
of the ABM, recording r a at each t £ {0,0.01,0.02,..., 1}, and then find the a that minimizes the difference 
(in L 2 norm) between the ABM realizations and the solution to Equation ^ (see Appendix [E]). The results 
of this procedure for ABM simulations performed on the physical contact and Facebook social networks with 
parameters r 0 = 0.25, 9 = 0.15, P = 0.3, c\ = 1, and C 2 = 9 are given in Fig. [5| For these parameters we 
find the optimal a is given by a = 0.835 and a = 0.834 for the Facebook and physical contact networks, 
respectively. Equivalently, the fitted value for the visibility threshold 1 — a is given by 1 — a = 0.165 
and 1 — a = 0.166, respectively. Holding the initial condition r 0 = 0.25 and parameters /3 = 0.3, Ci = 
1 and C 2 = 9 constant, we are able to repeat the fitting procedure to find fitted values of a for 9 £ 

4 For simulations of the ABM an initial condition r a (0) = rg is implemented by choosing pro .V] nodes uniformly at random 
to be active. 
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Figure 6: Parameter Region IIIO. Result of finding a as a function of 9 by fitting Equation ([3| (the SVF) 
to rep = 100 ABM simulations with parameters /? = 0.3, c\ = 1, C 2 = 9, and ?'o = 0.25 (Region IIIO, 
c* < 1 — at < 0). (Blue circles) Facebook social network. (Red diamonds) Physical contact network. Line of 
best fit for networks given by solid (m = 0.690 ± 0.022, b = 0.058 ± 0.008) and dashed (m = 0.764 ± 0.032, 
b = 0.050 ± 0.006) lines, respectively. 


{0.105,0.11, 0.115,..., 0.25}. These results are displayed in Fig. [6] and illustrate that for a large range of 
individual thresholds 9 the Facebook social network has a slightly lower fitted visibility threshold 1 — a than 
does the physical contact network. 

Observe that, for the parameters employed in Figs. 5][6 (9 = 0.15, /3 = 0.3, C\ = 1, and C 2 = 9), the 
SVF model lies in the stable police state dynamical regime (Region IIIO, c* = ci/(ci + C 2 ) < 1 — a < /3). 
We note that this is critical to the fitting procedure, since the parameter a appears in the solution to 
Equation © only when the parameters a , /?, ci, and C 2 lie in Region IIIO, see Appendix [E] We then also 
attempt to approximate the ABM by the SVF model in Regions Hie and III1 by using the fitted a displayed 
in Fig. [6j Figures |7flC)| show some representative results. As illustrated in Fig. [TJ for many combinations 
of parameters the SVF model is able to approximate the behaviour of the ABM in parameter Region Ille 
(left panels of Fig. [7]) and parameter Region III1 (right panels of Fig. [7|. However, we observe that under 
certain circumstances stochastic effects and details of the network structure become important and the SVF 
ODE model fails to approximate, even in a qualitative sense, the behaviour of the ABM. Specifically, the 
SVF model fails to capture the behaviour of the ABM when the initial condition is small (ro << 1), when 
the initial condition is near the boundary between the basins of attraction of r = 0 and r = 1, and when 
parameters ci, C 2 , a, and /3 are close to the boundary between parameter Regions IIIO, Ille, or III1 in Fig. 
|4a} see Figs. ii and |10| respectively. These are well-known limitations of population-level ODE models 
that attempt to approximate complex dynamics El- Note that these plots also provide comparisons with 
the binomial visibility function (BVF) model, which we introduce in Section |4j 

In the previous section we observed that setting v(r a ; 9) = v s (r a \a(9)) in Equation ([8]) exactly recovers the 
SVF model. Therefore, a natural strategy for improving on the ability of the SVF model to approximate 
the aggregate behaviour of the ABM is to specify an alternative visibility function that better captures the 
structure of the underlying network. Section [4] proposes two such visibility functions and compares their 
ability to approximate the ABM to the ability of the SVF model to approximate the ABM. 
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(a) Facebook social network Gf (1 — a = 0.131) 




(b) Physical contact network Gp (1 — a = 0.134) 


Figure 7: (Left) Parameter Region Hie (02 = 3, 1 — a < c* < j3). (Right) Parameter Region III1 (02 = 1, 
1 — a < f3 < c*). Time traces of rep = 100 ABM simulations (grey), solution to the SVF model (solid black), 
and BVF model (dashed black) with parameters 9 = 0.11, /3 = 0.3, Ci = 1, and tq = 0.2. 
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Time (t) 



(a) Facebook social network Gf (1 — a = 0.131) 



Time (t) 

(b) Physical contact network Gp (1 — a = 0.134) 



Figure 8: The SVF model can fail to approximate the ABM when (left) the initial condition is small 
(r 0 = 0.01), or (right) the initial condition is near the boundary between the basins of attraction of r = 0 
and r = 1 (ro = 0.1 w 1 - a). Time traces of rep = 100 ABM simulations (grey), solution to the SVF 
model (solid black), and solution to the BVF model (black dashed) with parameters 9 = 0.11, = 0.3 and 

Ci = C 2 = 1 (Region III1, 1 — a < (3 < c*). In the right panels, the BVF model performs much better than 
the SVF model. 
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Figure 9: The SVF model can fail to approximate the ABM when the initial condition is near the boundary 
between the basins of attraction of r = 0 and r = 1, i.e. when rg ~ 1 — a. (Left) Facebook social network 
Gf (1 — a = 0.179) and (right) physical contact network Gp (1 — a = 0.180): Equation ([3]) has solution 
in Region Hie, while ABM has solution in Region 1110. Time traces of rep = 100 ABM simulations (grey), 
solution to the SVF model (solid black), and solution to the BVF model (dashed black) with parameters 
9 = 0.17, (3 = 0.3, Ci = 1, and ci = 3, and initial condition r 0 = 0.2 (Region Hie, 1 — a < c* < 0). The BVF 
model performs much better than the SVF model. 




Figure 10: The SVF model can fail to approximate the ABM when parameters are close to the boundary 
between regions. (Left) Facebook social network Gf (1 — a = 0.163) with parameters in Region III1 (c 2 = 2.3, 
1 — a < < c*): Equation ([3]) has solution in Region III1, while ABM has solution in Region Hie. (Right) 

Physical contact network Gp (l — a = 0.165) with parameters in Region Hie (c 2 = 2.4, 1 — a < c* < (3): The 
SVF has solution in Region Hie, while ABM has solution in Region III1. Time traces of rep = 100 ABM 
simulations (grey), solution to the SVF model (solid black), and solution to the BVF model (dashed black) 
with parameters 9 = 0.15, j3 = 0.3, and C\ = 1 and with intial condition ro = 0.2. The BVF model performs 
much better than the SVF model. 
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4 Alternative Visibility Functions for Population-Level ODE 


Due to the difficulty of analyzing and computational cost of simulating the ABM, as well as the limitations 
of the SVF in approximating the ABM, we propose two alternative choices for visibility functions (different 
from the step-visibility function v s (r; a)) that can lend insight into how the ABM process behaves on different 
networks. Specifically, Sections |4.1|4.2| introduce the binomial and empirical visibility functions, respectively 
(BVF and EVF). These two visibility functions can be seen to be equivalent in the limit of large network and 
sample size, as shown in Appendix [F| Section 4.3 demonstrates numerically that the use of the BVF/EVF 
in Equation ([8]) represents an improvement over the use of the step visibility function v s (r;a) (SVF) when 
approximating the ABM for real networks. 


4.1 Binomial Visibility Function 

Consider the following derivation. Suppose that v £ V is a node in the network G = G(V, E ) with degree k, 
and that the fraction of nodes active in the revolution is r. Also, assume that the states of v’s neighbours 
are active with probability r independently so that the probability of v having j neighbours active in the 
revolution is 


r 3 (l- r ) k ~ j . 


The probability that the fraction of v’s neighbours exceeds the linear threhsold 6 is, therefore, 


k 

E 


r 3 { 1 - r) k ~ 3 = BinCDF(fc - \9k ]; jfc, 1 - r), 


(9) 


where BinCDF(x; n,p) is the cumulative distribution function for the binomial distribution with n trials and 
probability of success p evaluated at x. This quantity can be used to approximate the function v(r\ 0) (the 
expected fraction of the population that can see the revolution). In lj we considered k in Equation © 
to be the average degree over all nodes and argued that Equation © is a sigmoidal function that can be 
approximated by step-visibility function v s (r;a) for some appropriate visibility parameter a (see Appendix 
[B] for a brief summary of this argument). One possible improvement on adopting v{r\9) = v s (r;a{9)) in 
Equation ^ (resulting in the SVF) is to actually adopt v{r\6) = BinCDF([fcJ — \9k\, |_fcj, 1 — r) instead, 
where k is the average degree of the network. However, we can improve on this even further in terms of 
using more relevant information about the real graph by choosing v(r ; 9) to be equal to 


k— \0k~\ 


Vb 


JL1 1 //A 

(r;9,p k ) = ^2 p k BmCDF(k - \9k];k, 1 — r) = ^ E Pki .)0- ~ r) 3 r k ~ 3 . 

k k j =0 


( 10 ) 


We call the function Vb(r; 9 , p) in Equation (10) the binomial visibility function (BVF), see Fig. 11 Note the 
steep sigmoidal shape of Vb(r;9,pk) for the two empirical networks of Fig. El further justifying the choice 
of a step visibility function as an appropriate approximation in [l]. We note that computing the binomial 
visibility function on some grid of p + 1 discrete r-values r £ { r i}i=i requires calculating a double sum for 
every i = 1,...// + 1. In practice, this can be very expensive, especially for degree distributions pk with fat 
tails. For this reason, the following section introduces a complementary but less expensive approximation to 


the visibility function that is equivalent to the binomial visibility function of Equation (101 in the limit of 
large network and sample size. 


4.2 Empirical Visibility Function 

We now propose to reconstruct the binomial visibility function Vb(r;9, pk) empirically by sampling the 
underlying network directly. Specifically, for fixed 9 we approximate Vb by a discrete function v e £ 
where Vj = 1,..., p + 1 : v e j ~ Vb{^p-)- Since Ub(0) = 0 and Ub(l) = 1, we set v e p = 0 and r' e ,/i+i = 1- We 
calculate v e j for j = 2,..., p as follows. 
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(a) Facebook social network Gf (b) Physical Contact network Gp 

Figure 11: Binomial visibility function (solid green) and empirical visibility function (dashed blue) for Face- 
book and physical contact networks generated for linear thresholds 9 = {0.13,0.17,0.21,0.25} (curves from 
left to right). The binomial and empirical visibility functions are computed on the grid r £ {0, gjj, ..., 1}. 


1. Seed the network with r : j = -— active nodes in expectation. (Generate N random numbers drawn 
from the uniform distribution on [0,1]: randi for i = 1,..., N. If randi < then node i is active, 
otherwise node i is inactive.) 

2. Determine the fraction of nodes that can see the revolution. 

3. Repeat 1-2 a total of rep times and set v e j to be the average over the realizations. (We choose 
rep = 100 in our numerical experiments.) 


We denote the linear interpolation between the v e> j by the function i> e (r; 9) 


v e (r;9) 


J e,j 


J e,j 


j/l*-r 

1/m 


+ v eJ+l 


r-{j- 1)/V 

1/fj, 


if r = *^4 for some j = 1,..., fi + 1 
if r £ [^, i] for some j = 1, — , /x 


( 11 ) 


which we call the empirical visibility function , see Fig. EH 

We now make several brief remarks about the procedure for calculating the empirical visibility function, 
outlined above. Firstly, we remark that this procedure is easily extended to the case where 9 V are specified 
for each individual node. However, to keep our analysis as simple as possible we continue to restrict ourselves 
to the case where 9 V = 9. In Appendix [F] we show that the empirical and binomial visibility functions are 
equivalent in the limit of large network and sample size. Figure [TT] shows that they can also be expected to 
give similar results for finite network and sample size. Finally, we remark that compared to the binomial 
visibility function, the procedure for calculating the empirical visibility function is in practice significantly 
less costly to implement, especially for networks with fat-tailed degree distributions. 


4.3 Comparison of ABM and BVF/EVF Models 


In this section we compare the ability of the BVF model to approximate the ABM with that of the SVF 
model. We note that since the BVF and EVF are equivalent in the limit of large network and sample size 


and are very close in practice (see Fig. Ill, the observations we make about the behaviour of the BVF model 


also apply to the behaviour of the EVF model. We further note that the behaviour of the BVF/EVF model 
has parameters regimes that induce behaviour analogous to those observed in the SVF model. Since the 
BVF/EVF have sigmoidal shape, this follows directly from the analysis presented in Appendix B of (TJ for 
visibility and policing functions of general sigmoidal shape (the results from this analysis are summarized in 
Appendix [C]) . 
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We observe four types of outcomes when comparing the ability of the BVF model to approximate the ABM 
to the ability of the SVF model to approximate the ABM: 


1. The BVF model may succeed better in approximating the dynamics of the ABM than the SVF model, 
even if the improvement is only qualitative in nature, see Figs. [8^, and b (right panel), [9j and 10 (left 
panel), 

2. The BVF model may complement the dynamics predicted by the SVF model, see Fig. |12b[ 

3. The BVF model may produce qualitatively similar predictions to the SVF model, see Figs. 0 §> (left 
panel) and 10 (right panel), or 

4. Both the BVF and the SVF models fail to approximate the dynamics of the ABM either qualitatively 


or quantitatively, see Figs, (left panel) and 12a 


For each of these outcomes, the approximation of the ABM by the BVF model is no worse than the approx¬ 
imation by the SVF model, and, in many cases the approximation of the ABM by the BVF model is much 
better than the approximation by the SVF model. Therefore, we conclude that the BVF model, and hence 
the EVF model, represents an improvement over the SVF model with respect to their ability to approximate 
the behaviour of the ABM. This is no surprise, since the BVF/EVF model takes important information of 
the real network into account. 



(a) Both SVF and BVF models fail to approximate the 
dynamics of the ABM on the Facebook social network 
Gf (1 — &ode = 0.163) with parameters in Region 1110 
(c 2 = 1.4, 1 - a ode < 0 < e* « 0.417). 



Time ( t ) 

(b) The SVF and BVF models provide complementary 
approximations of the dynamics of the ABM on the phys¬ 
ical contact network Gp (1 — &ode = 0.165) with pa¬ 
rameters in Region 1110 (c 2 = 1.2, 1 — a ode < c* ~ 
0.455 < 0). 


Figure 12: The SVF and the BVF models can (a) both fail to approximate the dynamics of the ABM, or 
(b) act as complementary approximations to the dynamics of the ABM. Time traces of rep = 100 ABM 
simulations (grey), solution to the SVF model (black), and solution to the BVF model (dashed black), with 
parameters 9 = 0.15, /3 = 0.3, and Ci = 1 and with intial condition r$ = 0.1. 


5 Higher Order ODE Models 

In the previous section we introduced an extension to the SVF model. In particular, we showed that 
the BVF/EVF model is better able to approximate the ABM than the SVF model. However, both the 
BVF/EVF and SVF models are one-dimensional ordinary differential equation. This section now compares 
the performance of the BVF/EVF models as approximations to the ABM to the performance of a higher order 
(i.e., higher dimensional) compartmental model which we modify from [141 . We call this model the degree 
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approximation (DA) because it compartmentalizes individuals based on their state as well as their degree. 
We demonstrate numerically that the BVF model, and hence the EVF model, is no worse at approximating 
the aggregate behaviour of the ABM than the DA. This highlights the usefulness of the BVF/EVF model, 
especially since the DA is in practice much more computationally expensive to solve and much harder to 
analyze than the BVF/EVF model. 

5.1 Degree Approximation Model 

Before we detail the approach which we modify from El. it is useful to introduce some notation. We define 
14 to be the set of nodes with degree k and N k to be the set of nodes that have at least one neighbour of 
degree k, 


V k = {v GV : d v = k}, and 

N k = {w G V : 3 v GV k such that (u, w) £ E}, 

respectively, where we denote the degree distribution of N k by p k j, i.e. 

number of nodes of degree j that have at least one neighbour of degree k \Vj DN k \ 
number of nodes that have at least one neighbour of degree k |TVfc| 

As mentioned above, the degree approximation aggregates individuals according to their state and degree. 
Thus, conditioning on the state of the system at time tg, we define 

r (k) (t) = r (k) (t\t 0 ) = E M*)M 

1 k ' vGV k 

to be the fraction of nodes with degree k that are expected to be in state 1 at time t. Since nodes and edges 
are neither created nor destroyed, the fraction of nodes with degree k that are expected to be in state 0 at 
time t is given by 1 — r^ k \t). Analogously, we define 

r d (t) = r d (t\t 0 ) = l y fcl r(fc) W 

k 

to be the fraction of nodes expected to be in state 1 at time t. 

Applying the notation introduced above with the rules characterizing the dynamics of the process we find 
an equation analogous to Equation Q 

Ar<*)(i) = r«(f) - r<*>(t 0 ) = <? (fc) (^o) - d^ k \t\t 0 ), (12) 

where the expected growth and decay of the fraction of active nodes is modelled by the nonnegative growth 
and decay functions g^ k \t) = </ fc ^(t|to) and d^ k \t) = d( k \t\to), respectively. In order to close the degree 
approximation we need to approximate the quantities g^ k \t) and d^ k \t) in terms of r^ k \t), {9 v } veV , and /3. 

For notational convenience we first suppress the time argument, and consider the case where 8 V = 6. In 
order to approximate g^ k \t) or d^ k \t) in terms of r^ k \ 9 , and /3, we are now required to make an assumption 
about how the states of nodes sharing a common neighbour are correlated. Specifically, as for the BVF and 
EVF models, for any fixed v G V k we assume that the states of any two neighbours of v are independent. 
This assumption implies that the probability of v having exactly j active neighbours is (^)n/(l — n k ) k ~\ 
where n k = Y2i Pk,i r ^ is th e fraction of nodes in N k with state 1. So, the probability of v having at least 
\9k~\ active neighbours is 

( fc V J fc (l - n k ) k ~i = 1 - BinCDF([6»fc] - 1; k,n k ) = BinCDF (jfc - \9k] ; k, 1 - h fc ), 

It follows that 

g^ k \t)=c i (l-r (fe) (f 0 )) BinCDF(fc — [0*1; k, 1 - n*)Ai + o(Ai). (13) 
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As in Equation ([5]), we choose 


d [k) (t) =c 2 r {k \t 0 ) p(r d (t 0 );/3)At. + o{At). 


(14) 


Combining equations (12 )-|l4|) and dividing by At gives 
Ar® 


At 


= ci (l-r (fc) (t 0 )) BinCDF(A: — \9k ]; k, 1 - n k ) - c 2 r (k \t 0 ) p(r (fc) (t 0 ); /?) + o(l). 


Now, taking the limit as At —>• 0 yields 


= ci(l — r^)BinCDF(fc — \9k ]; fe, 1 — nk) — c 2 p(r^; /3). 


(15) 





(a) Facebook social network Gf- Parameters ci = 1, /3 = 0.3, and (Left) 0 = 0.19, c 2 = 3.0, ro = 0.2 (Middle) 
0 = 0.15, c 2 = 2.2, r 0 = 0.2 (Right) 9 = 0.11, c 2 = 1.0, r 0 = 0.1. 





(b) Physical contact network Gp. Parameters ci = 1, /3 as 0.3, and (Left) 9 = 0.19, c 2 = 3.0, ro = 0.2 (Middle) 
9 = 0.11, c 2 = 3.0, r 0 = 0.1 (Right) 9 = 0.11, c 2 = 1.0, r 0 = 0.1. 


Figure 13: The BVF model and the DA produce qualitatively similar predictions. From left to right ABM 
simulations die out (Region 1110), reach an intermediate equilibrium (Region Hie), or reach the fully realized 
revolution (Region III1). Time traces of rep = 100 ABM simulations (grey), solution to BVF model (dashed 
black) and DA (solid black). 


5.2 Comparison of ABM and Degree Approximation 


In Section |4.3| we showed that the approximation of the ABM by the BVF/EVF model is no worse, and 
usually much better, than the approximation of the ABM by the SVF model. In this section we show that the 
approximation of the ABM by the BVF/EVF model is no worse than the approximation of the ABM by the 
DA given in Equation (151. Specifically, we observe that the BVF and DA model produce similar qualitative 
predictions for most choices of parameters 9, /3, ci, and c 2 , see for example Fig. [l3j Also, although this is an 
atypical result, we present one set of parameters where the BVF model outperforms the DA on the Facebook 
social network Gf , and where the BVF model complements the DA on the physical contact network Gp, 
see Fig. [14] We observe that, for the specific parameters and networks used in the simulations of Figs. 
HfTT the DA does not significantly outperform the BVF model. It was confirmed in extensive additional 
simulations that this occurs generically for large parts of the parameter space of the models. This is worth 
emphasizing once more, since the BVF model is much easier to analyze and much less costly to solve than 
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the DA. This supports our proposition that the BVF (or EVF) model is a powerful tool for approximating 
and analyzing the ABM and a general technique for formulating tractable models for spreading processes on 
social networks that take real network characteristics into account. 



Time (t) 


(a) Facebook social network Gf with parameters ci = 1, 
C 2 = 3, 6 = 0.15, and /3 = 0.3, with initial condition 
ro = 0.2. The BVF model is a better approximation of 
ABM realizations than the DA. 
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(b) Physical contact network Gp with parameters ci = 
1, C 2 = 1.2, 9 = 0.15, and /? = 0.3, with initial condition 
ro = 0.1. The BVF model and the DA provide com¬ 
plementary information on how the ABM realizations 
evolve. 


Figure 14: Example of parameters for which (a) The BVF model outperforms the DA on the Facebook social 
network, and (b) the BVF model and the DA complement each other on the physical contact network. Time 
traces of rep = 100 ABM simulations (grey), solution to the BVF model (dashed black) and the DA (solid 
black). 


6 The Basic Reproduction Number ( Rq ) and the Initial Slope of 
the BVF/EVF 


In epidemiological models the basic reproduction number Rq is defined to be the number of secondary 
infections caused by a single infected individual introduced into a population that is entirely susceptible [9]. 
Note that in this context exposure to one infected individual is sufficient for a new infection. If R 0 < 1 then 
an infected individual is, on average, unable to replace himself and the outbreak terminates on its own in 
short order. On the other hand, if Rq > 1 then an infected individual is, on average, more than able to 
replace himself and the outbreak may spread into a full-blown epidemic. We define the basic reproduction 
number for our linear threshold model of political revolution analogously: the basic reproduction number is 
the average number of individuals that become active in the revolution due directly to the introduction of a 
single active individual into a population that is otherwise completely inactive. In this section, we show that 
the basic reproduction number in our linear threshold ABM is related to the initial slope of the BVF/EVF 
by the equation 

R o = -^v , b (q-,0,Pk), ( 16 ) 

Cl + C2_ 

—C* 

and discuss efficient ways to compute Rq. 

We begin by deriving an expression for Ro. Suppose that we activate individual v £ V in the population, 
and suppose that individual v has degree k. Further, suppose that this individual has V neighbours who can 
see the revolution (i.e., the linear threshold criterion is satisfied) once individual v is activated. Let w £ V be 
a neighbour of v’s that can see the revolution, and let t\ and r 2 be the first arrival time of Poisson processes 
with rates Ci and c 2 , respectively. The probability that individual w becomes active before individual v 
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becomes inactive is 


/»oo 

P (ti < t 2 ) = / P (ti < /c|r 2 = t) P (r 2 = £) dt 

Jo 

/*oo 

= / [1 — exp(—cif)]c 2 exp(— C 2 t)dt 

Jo 


= Cl = c «. 

Cl + C 2 

Since we only consider activations that result directly from ids activation (and not activations caused by a 
combination of the activation of v and the subsequent activation of ids neighbours), for our purposes the 
behaviour of each of ids neighbours is independent. Thus, the expected number of ids neighbours which 
become active is Vc*. 

Recall that the social network G = G(V ,, E ) has N nodes, M edges, degree distribution pk and secondary de¬ 
gree distribution Pk,j- For the single activated individual v with degree k the expected number of neighbours 
who can see the revolution is 

E[V] = kJ2Pk,jhje<i}- 

3 

Averaging over all possible degrees k now yields the basic reproduction number 


R 0 = c*E[V] = 


Cl 


Cl + c 2 


J2 k PkPkjhje<i}=c* kpkPkj• 


k,j 


k =1 j —1 


It is now possible to prove the relationship given in Equation (16). Writing 

L$ -1 J oo 
Cl 


it suffices to show that 


R °~ c +c k PkPk,j, 

j =1 k=l 


jpj — y ^ kpkPkj, 

k =1 


since 


v b (r; 0 , p k ) = ^2 P*3inCDF(fc - \6k ]; fc, 1 - r) 
k= 1 

oo k— \Bk~\ , . 

= IZ Pk[A(X-r) j r k -i 

k =1 j=0 
oo 

= J2p k 


k -1 


r + ... + 


k — \9k~\ 


(1-r) 


L0" 


fc -r efc i r r efe i 


^b(0; 6>, p fc ) = fcpfc = Y k Pk- 


k =1 


fc=l 


(17) 


We note that together with Equation (16), Equation gives an exact and cheap way to compute Rq- 

From the undirected social network G = G(V,E) form the directed social network G' = G(V,E'), where 
E' = {(v,w) : {iqic} £ E}. In this case, the number of edges emanating from all nodes of degree k is kpuN. 
So, the number of edges emanating from nodes of degree k and incident on nodes of degree j is kpkNpkj. 
Thus, the number of edges incident on nodes of degree j is 

OO 

NY kpkPkj. 

fc= 1 
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Equivalently, the number of edges incident on nodes of degree j is 


•^in¬ 


putting these two expressions together, we find 

OO 

jpj = 'y ] kpkPk,j, 

*:=i 


which completes the proof. 

Intuitively, this makes sense since the number of nodes that can see the revolution when one node is activated 
can also be approximated by 


Nv b 



=o 


vb (^;8,pk) ~v b (0 ;d,p k ) 


N 


v b (0 ;9,p k ). 


7 Application: Linear Threshold Propagation on Online versus 
Offline Social Networks 

In this section we investigate experimentally the differences in spreading behaviour that occur under the 
linear threshold model for the spread of political revolutions when applied to some empirical online and 
offline social networks, searching for some initial quantitative evidence that political revolutions may be 
facilitated by the network structure of online social networks of social media. 


7.1 Network Structure: Online versus Offline Social Networks 

It is often assumed that the connectivity of modern online social networks was an important factor in the 
spread of political revolutions in the past decade, e.g., in the Arab Spring revolutions of 2011 nansuraiisi, 
while traditional offline social networks (using in-person physical contact, or phone or mail interaction for safe 
communication) featured a different connectivity structure that was often severely restricted by censorship 
of the regime. Unfortunately, representative samples of the offline, traditional communication networks that 
were in existence prior to the adoption of new media technologies are unavailable at sufficient scale for 
countries affected by the Arab Spring, or indeed for any country. As opposed to the online networks of social 
media which are, by their nature, digitally stored and available, the offline social networks of pre-Internet 
societies have not been recorded at scale, simply because it was impractical in terms of cost and effort. This 
is a serious roadblock when investigating the effects of the structure of new media networks on the dynamics 
of political revolution, and if one wants to compare with propagation on pre-Internet social networks, it is 
necessary to identify proxy networks that are likely to be reasonable approximations to pre-Internet social 
networks, in terms of network structure. 

In this section, we use simulations of our ABM model on two empirical networks as a starting point to 
investigate differences in propagation properties that may arise between online and offline social networks 
within the linear threshold propagation model. We choose the small physical contact network Gp between 
individuals from [2U] as a representative for offline social networks, and we choose the Facebook network Gp 
as a representative for online social networks. 


7.2 ABM on Facebook Gf and Physical Contact Gp Networks 

We now briefly examine the differences in how the political revolution spreads on the Facebook social and 
physical contact networks via direct simulation of the ABM. In particular, choosing parameters 6 = 0.1, 
/3 = 0.3, ci = 1, and C 2 = 0.2 in the unstable police state region (Region III1, 1 — a < /? < c*), we simulate 
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Initial Condition (ro) 


(a) Facebook social network Gf 


(b) Physical contact network Gp 


Figure 15: Average final size rfi na i of rep = 100 realizations of ABM simulations on Facebook and physical 
contact networks with ro € [0,0.05], 0 = 0.1, /? = 0.3, Ci = 1, and = 0.2 (Region III1, 1 — a < ft < c*). 
Realizations are allowed to run until time tfi na i = 20. Basic reproduction numbers Rq for Facebook social 
and physical contact networks are Rq = 1.12 > 1 and Rq = 0.35 < 1, respectively. 


rep = 100 realizations of the ABM for each initial condition r 0 £ {0,0.0015, 0.003,..., 0.0495} from time 
t = 0 until time t = 20. The average final size (at t = 20) for each initial condition ( rf ina i\ro ) = (r a (20))| ro 
is recorded and displayed in Fig. [15] which demonstrates that for these parameters the political revolution 
described by our ABM propagates to a much greater extent on the Facebook social network, i.e. the online 
proxy network, than on the physical contact network, i.e. the offline proxy network. 

As an important illustration of the usefulness of the basic reproduction number Rq we defined in Sec. [6j we 
find that this is consistent with a difference in Rq for the two networks: the basic reproduction number is 
Rq = 1.12 (above the value of 1 which is expected to be required for propagation) for the Facebook network, 
and Rq = 0.35 < 1 for the physical contact network. 

These results indicate that the offline social network is less conducive to spreading the revolution in the 
ABM than the online social network: it has a smaller basic reproduction number, and in simulations a 
larger initial population of revolutionaries is required to spread the revolution. This provides some initial 
quantitative evidence that the spread of revolutions under a linear threshold process may occur more easily 
on modern online social networks than on traditional offline networks. While this is an interesting first 
observation, the next section discusses limitations and further investigations that are required to address 
this intriguing but complex question more comprehensively. 


8 Discussion and Conclusion 

In this paper we developed a linear threshold agent based model (ABM) to model the spread of a political 
revolution in a dictatorial regime. We showed that this model is consistent with previous simple step visibility 
function (SVF) ODE model developed in pQ. Using the relationship between these two models as a template 
we developed a hierarchy of models of varying complexity that approximate the behaviour of the ABM, see 
Table [l] Of all the models we have identified, we find that the BVF and EVF models (models of moderate 
complexity) offer the optimal combination of low computational complexity (cost), ease of analysis, and 
ability to approximate the behaviour of the ABM. Specifically, we find that for most parameters and initial 
conditions the BVF/EVF model is better able to approximate average ABM behaviour than the SVF. Also, 
for most networks the BVF/EVF model is much less costly to solve and much easier to analyze than the 
degree approximation from m- Importantly, the analysis of the simple ODE models in terms of stability 
of solutions for various parameter regimes directly gives insight in the qualitative dynamics of the linear 
threshold ABM for the spread of political revolutions. 
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We extended the concept of the basic reproduction number i? 0 from epidemiology [9] to the linear threshold 
ABM, we showed how it is related to the slope of the empirical or binomial visibility functions at r = 0, and 
we provided efficient ways to compute or estimate it. Analogously to epidemiological models, when Rq > 1 
we expect the political revolution to spread, and when Rq < 1 we expect it to die out. Thus, computing 
this quantity for a network can give an indication of how the ABM will behave on that network without the 
need to perform simulations, as we have demonstrated for empirical networks. 

The Facebook and physical contact networks we consider as case studies, i.e. the online and offline proxy 
networks, provide initial support to the hypothesis that the adoption of online social media may facilitate 
the spread of political revolutions by effectively changing the connectivity structure of the population in a 
way that makes linear threshold spreading more effective. Specifically, we find that for certain parameters 
the online proxy network is more susceptible to the linear threshold spreading process than the offline proxy 
network. Moreover, we find that the different behaviour of these two networks is consistent with the basic 
reproduction number calculated for these two networks. 

In addition to studying the ABM and its approximations on more and larger online and offline networks, 
much work remains to be done in the actual modelling of the political revolution process. For example, 
throughout this manuscript we have assumed that the linear threshold of individuals is constant for the 
population, that the underlying communication networks are static, that the graph is undirected, and that 
the nodes that are initially activated are chosen uniformly at random. Each of these assumptions represents 
a major simplification of reality that needs to be addressed to study further aspects of the spread of political 
revolutions on online social networks. For example, consider that in our model individuals estimate the 
current participation in the revolution by sampling their neighbours. Since individuals with greater sample 
size, i.e. larger degree, can form more accurate estimates of the current participation in the revolution, 
they should be willing to join the revolution at a lower linear threshold 6 than individuals who are more 
uncertain in their estimate of revolution size. We would therefore expect the linear threshold to vary from 
individual to individual as a function of their network degree. Or, consider that one of the principal strengths 
of new media is the ability of individuals to search for both content and like-minded individuals. Thus, we 
would expect that the underlying communication network should be changing on the same time scale as the 
revolution. Finally, one might consider that the nodes that are predisposed to be active at the initial stages 
of a revolution, e.g. nodes that represent activists, may also have larger degree (if, for example, they are 
charismatic and allowed to accrue followers) or smaller degree (if, for example, they are the specific target 
of regime censorship). There is significant value in the kind of parsimonious ABM model on static networks 
that we have considered in this paper, because much can be learned from this type of model and it is easier 
to analyze and interpret than more complex models. Nevertheless, extensions along the lines sketched above 
are important avenues for further study. 

The simulations in this paper were performed in Matlab. All Matlab code and data files will be made 
publicly available online upon publication of this paper. 


References 

[1] J.C. Lang and H. De Sterck. The Arab Spring: A Simple Compartmental Model for the Dynamics of a 
Revolution. Mathematical Social Sciences , 69:12-21, 2014. 

[2] D. Kempe, J. Kleinberg, and E. Tardos. Maximizing the spread of influence through a social network. 
In Proceedings of the ninth ACM SIGKDD international conference on Knowledge discovery and data 
mining , pages 137-147, 2003. 

[3] D. Centola, R. Wilier, and M. Macy. The Emperor’s Dilemma: A Computational Model of Self-Enforcing 
Norms. American Journal of Sociology, 110(4):1009-1040, 2005. 

[4] D. Centola, V.M. Eguiluz, and M.W. Macy. Cascade dynamics of complex propagation. Physica A, 
374:449-456, 2007. 

[5] D. Centola and M.W. Macy. Complex Contagions and the Weakness of Long Ties. American Journal 
of Sociology, 113(3):702-734, 2007. 


24 


[6] D.J. Watts. A simple model of global cascades on random networks. Proceedings of the National 
Academy of Sciences, 99(9):5766—5771, 2002. 

[7] G.M. Ames, D.B. George, C.P. Hampson, A.R. Kanarek, C.D. McBee, D.R. Lockwood, J.D. Achter, 
and C.T. Webb. Using network properties to predict disease dynamics on human contact networks. 
Proceedings of the Royal Society B: Biological Sciences, 278(1724):3544-3550, 2011. 

[8] S. Bansal, B.T. Grenfell, and L.A. Meyers. When individual behaviour matters: homogeneous and 
network models in epidemiology. Journal of The Royal Society Interface, 4(16):879-891, 2007. 

[9] H.W. Hethcote. The mathematics of infectious diseases. SIAM Review, 42(4):599-653, 2000. 

[10] P. Lind, L. da Silva, J. Andrade, and H. Herrmann. Spreading gossip in social networks. Physical 
Review E, 76(3):036117, 2007. 

[11] L. Zhao, Q. Wang, J. Cheng, Y. Chen, J. Wang, and W. Huang. Rumor spreading model with consid¬ 
eration of forgetting mechanism: A case of online blogging LiveJournal. Physica A, 390(13):2619-2625, 
2011 . 

[12] M. Granovetter. Threshold models of collective behavior. American Journal of Sociology, 83(6):1420- 
1443, 1978. 

[13] T. Kuran. Now Out of Never: The Element of Surprise in the East European Revolution of 1989. World 
Politics, 44(l):7-48, 1991. 

[14] M. Nekovee, Y. Moreno, G. Bianconi, and M. Marsili. Theory of rumour spreading in complex social 
networks. Physica A: Statistical Mechanics and its Applications, 374(l):457-470, January 2007. 

[15] Nahed Eltantawy and Julie B Wiest. The Arab Spring — social media in the Egyptian revolution: 
Reconsidering resource mobilization theory. International Journal of Communication, 5:18, 2011. 

[16] Habibul Haque Khondker. Role of the new media in the Arab Spring. Globalizations, 8(5):675-679, 
2011 . 

[17] Zeynep Tufekci and Christopher Wilson. Social media and the decision to participate in political protest: 
Observations from Tahrir Square. Journal of Communication, 62(2):363-379, 2012. 

[18] S. Khamis and K. Vaughn. Cyberactivism in the Egyptian Revolution: How Civic Engagement and 
Citizen Journalism Tilted the Balance. Arab Media and Society, 14, 2011. 

[19] Julian McAuley and Jure Leskovec. Learning to discover social circles in ego networks. In Advances in 
Neural Information Processing Systems 25, pages 548-556, 2012. 

[20] Marcel Salathe, Maria Kazandjieva, Jung Woo Lee, Philip Levis, Marcus W Feldman, and James H 
Jones. A high-resolution human contact network for infectious disease transmission. Proceedings of the 
National Academy of Sciences, 107(51):22020-22025, 2010. 

[21] Attila Giirsoy and Murat Atun. Neighbourhood preserving load balancing: a self-organizing approach. 
In Euro-Par 2000 Parallel Processing, pages 234-241. Springer, 2000. 

[22] I.Z. Kiss, D.M. Green, and R.R. Kao. The effect of contact heterogeneity and multiple routes of 
transmission on final epidemic size. Mathematical Biosciences, 203(1) :124-136, 2006. 

[23] R. Pastor-Satorras and A. Vespignani. Epidemic dynamics in finite size scale-free networks. Physical 
Review E, 65(3):035108, 2002. 

[24] J. Lindquist, J. Ma, P. Driessche, and F.H. Willcboordse. Effective degree network disease models. 
Journal of Mathematical Biology, 62(2):143-164, 2010. 


25 



A Physical Contact Network Data 


In this appendix, we briefly discuss the network extraction protocol for the physical contact network Gp 
from PHI] . 

The network presented in [20] was constructed by distributing wireless sensors to students, teachers, and 
staff at a U.S. high school during a one day period from approximately 08:00 to 16:30. When two wireless 
sensors are in proximity of one another, i.e. when they are less than approximately 3m apart, they register 
an interaction with a temporal resolution of 20s. Data are given for each seperate interaction in CSV format 
with three columns: ID of first wireless sensor, ID of second wireless sensor, duration of interaction (measured 
in 20s increments). A weighted undirected network Gpp is formed by connecting each pair of individuals 
with an edge whose weight is given by the total amount of time they spent in proximity to one another. In 
order to admit a comparison with the unweighted undirected Facebook subnetwork, we de-weight network 
Gp o by discarding all edges whose weight is less than the minimum duration wp, and by weighting all 
remaining edges equally. We denote the largest connected component of the resulting network by Gp, and 
for convenience we refer to this network as the physical contact network. 

The choice of the minimum duration wp is important since it determines, for example, the average degree of 
the physical contact network Gp. We choose wp such that the average degree of Gp is close to the average 
degree of Gf, for the following reason. In the simulations carried out in this paper to compare propagation 
under the linear threshold ABM on Gp and Gf , we seed the networks with a fixed percentage of active 
nodes and use the same threshold 9 on both networks. We want to calibrate the average degree of Gp to the 
average degree of Gp, such that, on average, nodes in Gp have the same chance as nodes in Gp to satisfy the 
linear threshold criterion and see the revolution. After this calibration, differences in propagation between 
Gp and Gf (for the same 9) are only due to the differences in network structure that go beyond the average 
degree. This motivates us to choose wp = 34 sensor measurements, equivalent to 11 minutes and 20 seconds, 
in order to match the average degrees of Gp and the Facebook subnetwork as closely as possible. 


B Justification of Step Visibility Function 


As argued in Appendix A of pQ, the expression for the visibility function given in Equation Q, with k 
taken to be the average degree of the network, can be expected to have a steep sigmoidal shape, which 


can be approximated by a step function. Indeed, Fig. 16 shows that BinCDF(LfcJ — \9k~\] |_fcj, 1 — r) has a 
clear sigmoidal shape, with a steep transition from 0 to 1. In |pQ it was also argued that increased average 
degree k can be expected to co-occur with decreased threshold 9 , since the more neighbors one has the more 
certain one can be about the true extent of the revolution because large samples are more reliable. Under 
this assumption, BinCDF(|_fcJ — ; |_fcj, 1 — r) can be considered a single-parameter sigmoidal function, 

which can be approximated by the single-parameter step visibility function (SVF) v s (r;a ), with parameter 
a indicating where the step transition occurs. 


We conclude this section with a few remarks regarding the assumption made in Sec. 4.1 and in [T] that the 
states of an individual’s neighbors are uncorrelated. While this assumption is applied in many models of 
spreading processes, for example rumour spreading processes [14] and epidemiological processes [22l 123] [24] . 
this assumption was not verified in [I]. In this paper we show that this assumption results in a valid 
approximation by comparing simulation results of the SVF model, and also the BVF/EVF models (which 
were derived in Sec. [4] under the same assumption), with ABM simulations for empirical networks (which do 
not make this assumption), finding consistent results for large parts of parameter space. 
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Figure 16: Dependence of BinCDF(LfcJ — |"0fc"|; |_fcj, 1 — r) on (top left) 0, (top right) k, and (bottom) 
negatively correlated 9 and k. 

C Analysis of BVF/EVF Model 


In this section, we summarize the parameter regime and solution stability analysis that was developed in 
Appendix B of [3] for the one-dimensional ODE model 


dr 

dt 


= f = C\ (1 — r) v(r) 
7(0 


c 2 r p(r), 



(18) 


with generic sigmoidal visibility function v(r) as in Fig. 17 Since the newly derived BVF/EVF in this paper 
have sigmoidal shape, the stability analysis of |1] applies to the BVF/EVF models. 



Figure 17: Visibility function v{r) with a fast transition from 0 to 1 that follows a sigmoidal profile. 


As is shown in [T], the dynamics of model ( |18| ) with generic sigmoidal visibility and policing functions 
closely follows the corresponding dynamics of model © in Region II (an open interval of equlibria C (0,1)), 
Regions IIIO and III1 (one equilibrium € (0,1)), and Region Hie (three equilibria £ (0,1)), thus establishing 
the equivalence in terms of dynamic behaviour of models (|3| (with step functions) and (18) (with sigmoidal 
functions, including the BVF/EVF models). We summarize the phase portraits of the different regions for 
model (18) in Fig. 18 which are analogous to the phase portraits of the SVF model shown in Fig. |4b| See 
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Appendix B of [T] for a detailed derivation and further explanation. For the BVF/EVF models, the values 
of equilibrium quantities like r* and r** depend on 9, /3 and the network structure. 
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Figure 18: Equilibria, stability and basins of attraction on the r-axis (r € [0,1]) for model (18). Closed (open) 
circles represent locally asymptotically stable (unstable) equilibria. Left (right) arrows indicate regions where 
f < 0 (f > 0). The phase portraits are similar to the case of the SVF model in Fig. |4b| 


It is also important to point out how the SVF stability region plot of Fig. |4a| generalizes to the BVF/EVF 
models. For a given empirical network, the BVF/EVF model does not feature a range of visibility parameters 
a like the SVF model, but each choice of 9 results in a different BVF/EVF function. For a fixed 9 , the 
different phase portrait behaviours of Fig. [l 8 ]occur depending on the values of f3, C\ and c 2 . As a result, for 
a given empirical network, the parameter space for the BFV/EVF models as a function of 9 , /?, c\ and c 2 can 
be divided in regions as for the SVF model in Fig. |4a[ where 9 replaces a as the parameter on the horizontal 
axis. The precise shape of the BVF/EVF parameter regions would depend on the actual empirical network 
considered. For example, the boundary between the II and III regions is in general not a straight line, but 
needs to be determined numerically. 


D Gillespie’s Algorithm 

Following [7j, we numerically simulate the ABM by implementing Gillespie’s Algorithm. In this appendix 
we give a brief overview of this algorithm. We first introduce some notation. As above, let r a denote the 
fraction of individuals in the population that are expected to be active in the revolution at time t. For each 
node v £ V we let j v (t) = 1 if v can see the revolution at time t, i.e. if v satisfies Equation ([I]), and let 
7 v (t) = 0 otherwise. Furthermore, we let £-| ,v and £ 2 ,v denote the first arrival times of independent Poisson 
processes with rates c\ and c 2 , respectively. 

Gillespie’s Algorithm is based on the fact that the sum of two independent Poisson variables is also a Poisson 
variable with rate equal to the sum of the rates of the original processes. It follows that 

Y Km “ Sv (*)) TuW + s v (t) p(r a {t )\/?)] 

v£V 


is a Poisson process with rate 


A = Y, [ci (1 ~ s„(f)) 7 v (t) + c 2 s v (t ) p(r a (t)\P)] ■ 

vEV 

The first arrival time of this process, therefore, is an exponential random variable r with rate A. At time 
t + t the state of exactly one of the nodes will change. Moreover, since and £ 2) „ are independent, the 
probability that the state of node v will change is 

ro G (1 s v (t)) 7 v (t) + c 2 s v (t) p(r a (t);/3) 

A 

The Gillespie Algorithm then proceeds iteratively in three steps. 
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1. Find the time r of the next event by drawing r from an exponential distribution with rate A. 

2. Determine which node changes state by drawing one node from V = where node v € V is 

drawn with probability P^. 

3. Update t <— t + r and re-calculate s v (t), r(t), 7 „(f). 


E Complete Analytic Solution to SVF Model 


Since the complete analytic solution to the SVF model was not provided in [T] but is useful in the current 
paper, we provide it here. 

The solution to Equation ©> is as follows. In Region I (1 — a = (3) and Region II (1 — a > j3) 

( roe _c A t_t °) if r 0 < p 

r{t) = < r 0 if i"o £ [P,1- a) . 

[ 1 — (1 — ro)e -Cl ( t- *°) if r 0 > 1 — a 

In Regions IIIO (c* < 1 — a < f$). Hie (1 — a < c* < (3), and III1 (1 — a < f3 < c*) 


where 


r 0 e 


-c 2 (i-io) 


if rg < 1 — a 


[c* + (r 0 - c*)e-^+ c ^ t - t ^] I {t<ta} 

+(1 — if r 0 S [/?, 1 — a] and c* < 1 — a < /3 


r(t) = c* + (r 0 - c *) e ^( Cl + C2 )(*- to ) 


if rg € \/3 ,1 — a} and 1 — a < c* < /3 


[c* + (r 0 - c*)e-( Cl+C2 )(‘- t0 >] I {t<tp} 

+ [l - (1 - /3)e _Cl(t_t ' 3) ] if e [P, 1 - a] and 1 - a < < c* 


1 - (1 

toe — to 

tp = to — 


1 


Cl + C2 

1 

Cl + C 2 


log 

log 


if r 0 > /3 


1 - a - c* 
r 0 -c* 
p-c* 
r 0 - c’ 


, and 


F Equivalence of Binomial and Empirical Visibility Functions 


We briefly argue that the empirical and binomial visibility functions are equivalent in the limit of large 
network and sample size, i.e. as N,rep oo, see also Fig. EH Let v £ V be a node in the network 
G = G(V,E ) with degree k. Now, fix 9 and suppose that we are calculating v e j ~ = Vb{rj) via the 

algorithm presented in Sec. 4.2 The I th iteration of Step 1 of this algorithm selects A/j nodes uniformly at 


random to be active and results in Vi nodes that can see the revolution. Note that by the Law of Large 
Numbers 


rep 

is 


N 

V; 


->• n = 


j -1 


as N —> oo, and 


rep N 


E 


= E 


W 

N 


as rep —>• oo, 
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where Y denotes convergence in probability. Since the A/j activated nodes are chosen uniformly at random, 
the statuses of the neighbours of v are independent. Therefore, the probability that v can see the revolution 
is 


E 

l— |" 6k~\ 


MY 


i \N 


M 


N 


k-j 


Y l-Y = BinCDF k — \dk ]; fc, 1 


M 


N 


It follows that the expected fraction of nodes that have degree k and can see the revolution is 

Mr 


E 


Pk BinCDF k — \8k] ;k, 1 — 


N 


and hence, the expected fraction of nodes that can see the revolution is 


E 


= E 


L k 


M 


E PfcBinCDF ( k - \9k ]; k, 1 - -Y 


We conclude the proof by observing that by the Continuous Mapping Theorem Y rj and N —> oo implies 


E 


and hence, 


N 

5>BinCDF (k - \9k];k, 1 - rj) = v b (r;9,p k ) as N —> oo, 

k 

v e j —> Vb{rj\ 6 , pk) as N —> oo, rep —> oo. 
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